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Abstract : We present a strong convergence implicit Runge-Kuttci method, with four stages, for solution of 
initial value problem of ordinary differential equations. Collocation method is used to derive a continuous 
scheme; and the continuous scheme evaluated at special points, the Gaussian points of fourth degree Legendre 
polynomial, gives us four function evaluations and the Runge-Kutta method for the iteration of the solutions. 
Convergent properties of the method are discussed. Experimental problems used to check the quality of the 
scheme show that the method is highly efficient, A - stable, has simple structure, converges to exact solution 
faster and better than some existing popular methods cited in this paper. 

Keywords: strong convergent, collocation method, Gaussian points, highly efficient, A - stable, simple 
structure, 

I. Introduction 

A number of practical and real-life problems are models of ordinary differential equations (ODEs). Most of 
these are non-linear, stiff or oscillatory problems. Examples of such are found in Electric circuits systems, 
chemical reaction problems, population growth, Biology, Social Sciences etc. Explicit Runge-Kutta methods 
cannot efficiently handle them as they lack sufficient stability region. The instability of explicit Runge-Kutta 
(RK) methods motives the development of implicit methods (see 5). The coefficient of explicit method is a 
lower triangular matrix and that of the implicit is a square matrix. We have different types of implicit methods; 
Radau, Lobatto, Gauss Legendre quadrature. We have the Gauss-quadrature method of order four and six 
respectively (see 2). In this paper, in other to get higher and more stable method, we used four Gaussian points 
of Legendre polynomial to derive a four stage Gauss quadrature R-K method of order 8.which more efficient 
and easy in their implementation than the existing ones. 

" ri l'- 

II. Basic Definitions/Preliminaries: 

(i) A Runge-Kutta implicit scheme for first order ODEs is defineds as 

5 

Y n +1 = y n + h Fj, i = 1,2,... S ( S = 

i=l 

5 

where = f(x n + c ih , y* ), y„ + ^ ^ fj 

i=l 

(ii) Order of scheme: the order of implicit Runge-Kutta methods whose abscissae or integration Rodes 
are Gaussian points of Legendre polynomial of order p = ( 2S) 

(iii) Error/error construct: 

A Runge-Kutta solution can be expanded into Taylors series, the solution agrees with the Taylors 
series expander up to order P, the truncation error is the O(h p+1 ) term. The error and error constant 
is associated with a linear difference operator L(y(x), h). as L(y(x), h) = y(x n+1 ) — y n+1 


stages) 


/ i— ", 

► 

1.01 
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= Co y(x n ) + Cihyi (x n )+ (2h 2 y n (x n ) + ... C p hPy ( P } (x n ) + C p+1 h (n+1) y(x n ) 

where y(x n+1 ) and y n+1 are exact and Runge — Kutta approximate solutions, 

respectively. 

According to (5) the method is of order P, with error constant C p+1 of = C, = C 0 = C 2 C p = 
0 and C p+1 

(iv) The three stage implicit Gauss-quadrature (3) of order 6 is summarized by the table below: 


1A. Buchers Coefficient table (1) 



Gen. Solution 

5 

y n +1 = y n + h ^ b; f 

i—1 

5 

F = fix + c ih .yd.yi = y n + h ^ a i7 f l 

7=1 

See (8) 

III. Methodology 

We consider the font order initial value problem of ordinary differential equation. 
y 1 = f(x,y),y (x 0 ) = y 0 , a<x <b. - (2.01) 

We use a polynomial of the form 

t-l m -1 

y(x) = ^ (x)y n + h ^ Bj ( x)f{cj ) - (2.02) 

7=0 7=1 

where t are the number of interpolation points x n+ ,- (_/ = 0,1,... t — 1) m is the number of collocation 
points Cj () = 0,1, ...m — 1 ),/(x,y) is continuously differentiable. The constants Ky, Bj are elements of 
(t+m) x (t+m) square matrix. The are selected so that high accurate approximation of the solution of (2.01) 
is obtained, h is a constant step size. 
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The ( x),Bj (x)in (2.02) can be represented by polynomial of the form. 

t -1 

IX; (x) = ^K;,i +1 x‘ , j = (0, 1, ... t - 1) 


i —0 


2.03 


(x) = ^ h Bj i+1 x 1 , j = (0,1, ...m — 1) 

7=1 

le coefficient Bj,i + 1 are to be detei 

m+t-l t-1 m—1 

y (x) = ^ (^!X;, i+ i y„+i + h ^ h Bj i+1 f n+] ) X ; = ^ X; x ; - (2.04) 


7=1 

The coefficient K; i+1 Bj, i + 1 are to be determined. Now substituting (2.03) into (2.02), we have 

m+t-l t-1 m—1 m+t-l 


7=0 7=0 


7=0 


Where 


t-i 


m—1 


>y = ( X/ K/ ' i+1 yn+ i + h X f n+ i )’ x 7 G Ri 'j G (0, 1 ... t + m — 1). 


.7=0 


7=0 


(2.04) can be exposed in matrix form as y(x) = (y n ,y n+ i,... y„+t-i,/„,/„+i - /n+m-i) C r D 
Where 


K 0 1 ... . 

■ K t-1, 1 

hB 0 1,... 

hB m - 1,1 

A 

ix 0 2 .... 

■ “t-i, 2 

hB 0 2,... 

hB m _ i,2 


IX 0 3 ... . 

■ K t-l, 2 

hB 0 3,... 

^m-1,3 


K 0,t+m 

■Xt-1, 

t+m hB 0i 

t + m,.... hP m _ 1 t + m 



V 


(2. 05) 


1 x n x„ 


D: 


-t+m—1 




0 1 2 x n+i?1 . (t + m - 1) 2 


+qi 
t+m— 2 


0 1 2 x n+qm _ t ... (t + m- 1) x^J, 


(2.06) 




< 7 ; are the collocation points. 

Theorem 1.0 

Let I denote identity matrix of dimension (t + m), matrices C and D are defined by (2.05) and (2.06) 
respectively, satify. 

(i) DC = I 

(ii) y(x) = *' - (2.07) 

Now we assume a power series solution of degree order 4 of the form. 


y(x) = ^x 7 x 1 , y‘(x) = ^x ; t x 1 1 (2.08) 

7=0 7=0 

Interpolate at x n and collocate at x = x n+qi , i = 0,... 4, yields system of simultaneous equation of the form. 
X 0 + X; X n + x 2 xl + X 3 xl + \ 4 x% = y n 
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0 Xi+ 2 \ 2 X n+qi + .... 4 X 4 Xn +qi = fn+Rl 

0 X 4 2 X 2 X n+q 4 + 3 X 3 4 X 4 X n+q4 ~ fn Q 4 
Where X ; are to be determined. 

(2.09) can be rewritten in matrix form as 


r 1 

x n 

x ^ 

A/ n 

3 

x n 

0 

1 

9 Y 

3r 2 

0 

1 

9 r 

Zj - /t n+q , 2 

3r 2 

' JA -n+q2 

0 

1 

9 y 

■^^ 72+^3 

3r 2 

0 

1 

Y 

3r 2 

'J^-n+q 4 


V 

i.e 


v 4 X 



4-^n+qi 


X 1 

4*n+(?2 


X 2 

4^ n +q 3 


X 3 

4^ri+q4 


X 4 

J 


l J 


(2.09) 


r \ 
y n 

fn + <7l 

fn + <72 

fn + <73 


fn + <?4 


( 2 . 10 ) 


DA = Yj where 


A 



1 

x„ 

A, n 

3 

x n 

X 4 

A, n 

D= 

0 

1 

Y 

3r 2 



0 1 T V *3 -y 2 A y 3 

yu j. -^4 JA n+q 4 “ A n+q 4 A 


A = 


(x, 




Xq X 4 T 


D 


A — y n ,i 


^ fn+qi> 


fn+q ? i fn 


+93' 



( 2 . 11 ) 


Using maple mathematical software we obtain the continuous scheme of the form. 


i.e y(x) y n + h ( f n +q 1 + fn+q2 + fn+qj, + fn+q 4) 

y(x) = y n +r - <7 i<72<74* (h 2 q 2 q 4 + h 2 q 4 q 3 + h 2 q x q 3 + 

h 2 q^qf)x 2 _ _ 

(-94 + <7i)(-< 72 + <7i)(<7s + <7i) + 2 /i 3 C —<73 + <7i)(-< 72 + <7i)(<7 3 + <7i) 


- (q ? h + q 3 h + q 4 h)x 3 _ 

3 h 3 (—q 3 + q t )(-q 2 + <7i)(-<7s + <7i) 

fn+ < 7 i 


+ 


<7i<73<74* _ 

(<72-<?4)(-<72 + <7l)(<72 


73 ) 


+ 47t 4 (-q 4 + q 1 )(~q 2 + <7i)(—q 3 + <7i) 

(h 2 3 q 4 + h 2 ql q 4 + h 2 ql q 3 )x 2 
2h 6 {q 2 - q 3 ){~q 2 + qf)(q 2 - qf) 
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+ ( q 1 h + q 4 h + q 3 h)x 3 _ 

3 h 3 (q 4 - q 3 X~q 2 + qi)(q 2 ~ 
fn +q 2 


4 h 3 (q 2 - q 4 )(-q 2 + q 4 )(q 2 ~ q 3 ) 


+ 


q 4 q 3 q 4 x 


(q 3 - q 4 )(.q A q 2 ~ q x q 3 + q 3 ~ qiq 3 ) 

— (q 7 h + q 4 h + q-i h)x 3 _ 

3 h 3 (q 3 - q 4 )(q t q 2 - q x q 3 + q 3 ~ 
q2q3 fn+q3 


(h 2 al q 4 + h 2 al q 2 + h 2 2 q 4 ) x 2 
+ 2h 3 (q 1 q 2 -q 1 q 3 +q%-q 2 q 3 )(q 3 -q 4 ) 

_ x 4 _~" 

q 2 q 3 ) + 4 h 3 (q 3 - q 4 )(q t q 2 - q 4 q 3 + q 3 ~ 


+ 


qiq 2 q 4* ______ 

qiqiq-i - q\q 2 q\ + qfqi - q^q* + qiq 2 - q 2 q^q\ - ql+qiq 3 - 


(h 2 q,q ? + h 2 q,q 2 + h 2 q 7 qf)x 2 


2h 3 (q t q 2 q 3 - q 4 q 2 q 4 + q|<h - q 4 q 3 q 4 + qlq 2 ~ q 2 q 3 q 4 ~ ql + q%q 3 ) + 


(q-ih + q 7 h + 


3 h 3 (q 4 q 2 q 3 - q 4 q 2 q 4 + q\q 4 - q x q 3 q 4 + q\q 2 - q 2 q 3 q 4 ~ ql + qlqf) 
x 4 _ 

4h 3 (<M 2 93 - qiq 2 q4 + qlqi - qiqsq 4 + q\q 2 - q^q* -ql + qiq i) fn + Qa 

( 2 . 12 ) 

(l 3V206N 33 67 (1 , 3V206 \ 

Now evaluating the continuous at m =-, q 2 — —, q 3 — — ,q 4 — - H-, 

& V2 100 / 100 ^ 100 V2 100 ) 

We obtain discrete schemes: 


f 1633 71V206 \ 

y n+qi - y n + ( 18780 “ 96717000 


hfn- 


+ 


/134689 927a/206\ 


+ 


+ 91^^939000 78250 ) h ^ 

( 1633 121979V206\ 


n+q 2 


(171511 927a/ 206^ ^ f _ __, ^ 

\939000 ~~ 78250 ) ^ n+(?3 + \18780 19343400 J ^ 


n+q 4 


2n+q 2 Vn 3 ~ 


/ 7623 1629507a/206\ 


+ 


\78250 ' 257912000 >'" jn+q 1 

/ 3470313 118701 \ 


) 


hfn 


+ 


121284000 hfn+cn 4256800^' Vn+(?3 
/ 7623 1629507V206\ 




+ 




78250 257912000 


- hf r 


n+q 4 


y n + q 3 yn+ \ 117375 + 


/ 8978 1629507a/206\ 


257912000 


1 hf r 


n+q i 


+ 


4520423 

12770400 


hf 


n+q 2 


+ 


10410661 

63852000 


hfn+q 2 + 


( 8978 1629507V206^ 


hfn 


( 1633 

y n +q 4 - y n + I-+ 


+q3 \117375 257912000 1 Jn+q4 

121979 V206 \ /134689 927 V206 N 


\ 18780 ‘ 19343400 /^ n+(?1 + l 939000^ 78250 > "' Jn+q2 


+ 


hfn 
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+ 


/ 171511 
V939000 


927V206 

78250 


)hf. 


n +Q3 (l8780 


1633 17V206 


)hf. 


n+g 4 


- (2.13) 


96717000 . 

To convert to Runge-Kutta function evaluations, the discrete schemes ( 2.13) must satify 
(2.01), hence 

7 , . / 1633 71V206 \ 

y-n+qi ~ f\. x n+q 1 > Vn+qi) ~ f ( x n+q 1 >yn + I j 


+ 


/134689 927a/206\ 


y939000 78250 

/ 1633 
\18780 

y-n+q 2 ~ f(. X n+q 2 ’yn+q 2 ') 


jhfn+q 2 + ( 


18780 96717000/ -^ n+qi 

(171511 927a/206\ 


+ 


^939000 78250 

121979V206\ 

19343400 J h f n+ ^ 


- hf r 


n+q 3 


f (, x n+q 2 ’ Vn 3” I 


118701 


hf 


n+q 2 


+ 


/ 7623 




78250 


4256800 

yn+q 2 — f{, x n+q 2 ’yn+q 2 ) — 

f( x v + ( 8978 + 1629507 V206 N 

7 1 n+q 3 >yn 1,117375 257912000 / ' LJn 

10410661 . . ( 8978 


f 7623 i 1629507V206\ 
78250 + 257912000 

1629507a/206\ 
257912000 J ^ n+( ? 4 ] 






n+?l 


+ 


3470313 

21284000 


hf n 


+q 2 


+ 


hfr 


4520423 , . 

+£?1 + 12770400 


+?2 


+ 


63852000 '~ jn+q3 

Vn+q 4 — /(-^n+q 4 < yn+q 4 ) 


f (, x n+qAi Vn 3” 


\117375 

/ 1633 


1629507a/ 206\ 
257912000 j k ^ n+q4 


+ 


121979V206\ 


+ 


/171511 927a/206\ 


\939000 


+ 


78250 


7 




V18780 ■ 19343400 J 

/ 1633 




n+?i 


+ 


/927a/206\ 


\ 78250 


-) hf t 


n+q 2 


n+q 3 


+ 


+ 


71V206 ^ t ^ 
96717000 J hfn+q* 


yl8780 

Putting y n+(?1 — 7^1 — fi^n+qr’yn+qi)®^ fn+qi’K 2 ~ f{ x n+q 2 ’yn+q 2 )^ >x fn+q 2 6tC., 
We obtain the function evaluations as 


Ki =f\ x n+[ 2 


3a/206\ 

-Jh, f n + 


f 1633 71V206 \ 


100 


V 


18780 96717000/ 


hki 


+ 


/134689 927a/206\ 

1 -jhk 2 


\939000 


+ 


/171511 927a/206\ 

- hk 2 + 

1939000 78250 / 3 


78250 
/ 1633 

li 


121979a/206\ 
19343400 / hk ^ 


18780 

, 1629507a/206\ / 3470313 \ 

K 2 - f[ x n + Yo 0 h ’ yn + 1 78250 + 257912000 ) hk i + h 121284000/ hkl 


33 


/ 7623 


+ 


( 7623 

178250 


/ 118701 \ 

~~ V425680oJ hk3 

1629507V206\ 

- /i/c 4 ] 


257912000 


jr _ rr ,67 , / 8978 

K 2 — / fXv. H- h, y„ + (-h 

j J L n 4 qq ’ yn \ 11757c 


V117375 


1629507 V206' 
257912000 


\ , , , ( 4520423 \ , , 10410661 , , 

n/ci + - )hk 2 H- hk 2 

J 1 V12770400 / z 63852000 J 


(2.14) 


(2.15) 


Manuscript id. 371428231 


www.iistre.com 


Page 22 











































































Four Point Gauss Quadrature Runge - Kuta Method Of Order 8 For Ordinary Differential.. 


+ 




( 8978 1629507V206\ 

117375 257912000 

(1 3a/206\ 

*4=/k + 2+^or r y ” 


) 


hk 4 ] 


+ 


/ 1633 121979a/206\ 

l18780 + 19343400 / 


hki 


/134689 927a/206\ 

+ l 939000 + 78250 J hI<2 


+ 


/171511 927V206\ 




939000 


+ 


78250 


hk 3 + I 


/ 1633 71V206 \ 

18780 + 96717000 


)hk 4 ] 


The weight is b = (b t , b 2 , b 3 , b 4 ) ; Evaluating the continuous scheme (2.12) at x 
1 we obtain 


b = 


/ 1633 1531 1531 1633 N 

V9390 ’ 4695' 4695’ 9390, 


(2.16) 


The four stage Runge-Kutta method is defined as 

4 

1633 


Z 1633 „ 1531 „ 

bi ki = y n (/Ci + K 4 ) + ( K 2 + K 3 ) 

t =i 

Where K t i = 1,2, ...3 are given by (2.15) 

3.00 Analysis of the scheme. 

We can be write (2.15) in Butchers tableau as: 

A = (a y ), i,j = 1,2, ...4. 


1 

3V206 

( 1633 

71V206 ^ 

(134689 927V206\ /171511 

927V206\ ( 1633 

121979 V206^ 

2 

100 

V18780 

96717000 / 

V939000 78250 ) V939000 

78250 / V18780 

19343400 / 


33 

100 


67 

100 


/ 7623 _ 1629507 V206\ 347013 

V78250 257912000 ) 21284000 


118701 ( 7623 1629507 V206^ 


8978 1629507V206^ 


1 3V206 

2 ' 00 


(ll7375 

(■ 


257912000 


1633 121979 V206^ 


V18780 19343400 


4520423 

12770400 

/134689 
V939000 


4256800 

10410661 
63852000 Vll7375 


V78250 
(■ 


257912000 


( 2 . 16 ) 


8978 1629507 V206" 

257912000 


) 


927V206\ /171511 927V206\ / 1633 71V206 \ 

78250 ) V939000 + 78250 ) V18780 + 96717000 ) 


1633 

9390 


1531 

4695 


1531 

4695 


1633 

9390 


Or can be summarized as 


c 

A 

U 


B 

V 


Where (C = C X C 2 ... C 4 )f A = (a y ) y = 1,2, ...4, U = (1,1,1,1)f V = (1) 
(i) Consistency: 

Te Runge-Kutta method is consistent since 
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(ii) 



CLjj (* i, 



( see Table 2.16) 


Stability: 

The stability of the method is investigated by considering the linear test equation. 
y / = x y, x G <£ 

PuttingZ = xh, h G (0,1) 

The stability function is R (Z) 


R (Z) = I + Zb T (I — ZA) —1 


where I is the identity matrix 
b = b t> 

b 2 , b 3 , bf),the weight, e — (1,1,1,1), A is the Runge — Kutta matrix of 
coefficients (2.16). 

The region of A- stability is the set of point satisfying R (Z) = { Z: R (Z) 

< |) Ipid Re Z < 1} 

The characteristic polynomial is defined as 

P (Z) = det ( R(Z) -1 x) , which can be potted using maple or math lap, to 
determine the region of A- stability. 

(iii) Order and error of the scheme: 

The order of the method is p = 2S = 8, since the nodes or absyssaes are Gaussian 
points of fourth order Gauss - Legendre Polynomial and all Gauss- quadrantine 
methods have order P = 2S, ( S = stages) [3]. The Runge - Kutta solution y n+1 
can be expanded into Taylors series as 

i=1 i =1 
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4 

= Vn + h y "+«i 

i =1 

= Jn + C\hyL + C 2 h 2 yl l / + ... C hPy n{p+l) 

The R.K solution agrees with taylor’s series expansion up to term in h 8 . The truncation error is 
0(7i 9 ). 

Numerical Experiments / Results: 

To show the efficiency of this new method, we use three problems to compare our computational 
solutions with exact and other similar methods. 

Problem I. (Electric circuit problem model) 

In an Electric RL - circuit system, L is the inductance constant, I the current, R the resistance 
constant, and E9t), periodic electronmotive force. Find the currents at time t = 0.2 (0.2) 0.1 
secondly; with L = .1, R 5, E(t) = sin(t), with initial condition I (0) = 1.0. 

Model : By Kirchlotts voltage law (KVL) the voltage across the resister is RI, the voltage drop 
across inductor is E L . By KVL Law the sum of two voltage must be equal to electromotive force 
E(t). 

Thus our problem reduces to differential equation L ^ + RI = E(t) 
i.e. Y t = -50 1 + 1° sin(f), y(0) = 1.0. 

Example 2: (Logistic Population Model) 

A logistic mathematic model has useful application to human population and animal populations [ ]. 

The mathematical model is defined as 

y' = Ay{t) - B y 2 (t), y(t 0 ) = y 0 i A, B >0. 

A A A 

If y(0) < the population is monotonic increasing to limit - and if y(0) > it decreases to limit 
B B B 

A 

B' 

Now we consider a logistic problem 

y = Ay(t) — B y 2 (t), y(t 0 ) = yo with A = 5 6 = 3 initial condition 
y(0) = 2 Million Cows in Kaduna State, Nigeria, t is measured in months. 

Model: The logistic problem reduces to differential equation 
y = 5y — 3y 2 , y(0) =2, h = .2 

Example 3: (Chemical Reaction Model): 

A bimolecular reactions a and B combines a£ A moles per litre of distance A and b 6 6 moles of 
substance B given y(t) is the number of moles per litre that have reacted after time t we can find the 
number of moles per litre which reacted at given times. 
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Model: We use chemical reaction Law) which states that under constant temperature, the chemical 
reaction is proportional to the product of concentrations of the reacting substances, thus the model can 
be defined by a differential equation: 

y = K(a — y)(b — y), y(t 0 ) = Yo> K is the proportional constant, a, b are constants. 

In our example we consider a = 5, b = 3, y(0) = 0 and K = 1. Thus our model reduce to an initial 
value problem y = y 2 — 8y + 15, y(0) = 0. 

We can compute the solutions at t = 0 (.1) .5. 

The following are the solution tables for comparison. 

Notations: 

y(x n+1 ), n = 0,1,... = Exact solutions, with step — size h. 
y n+1 = The new R-.K approximate solutions. 

Abs = /y( n+ i) — y n+ i/, absolute error. 

ERi = absolute error of — point Gauss — quadrature (3)method. 

ER 2 = absolute error of Lie et al (8)method 


Table 1: Comparison of numerical solution of problem I. (Rl-Circuit). 


A 1 * 1 * Z' \ 500 ’ Z' \ , 2511 —KOt t0 . 

Analytic solution, y(t) =-sin(t) H- e -cosiTf) 

3 3 v 3 2501 v 3 2501 2501 v J 


t 

y(0 

y n +1 

ER 1 

er 2 

Abs 

.02 

0.3693509037 

0.3693509135 

3.81 E-06 

3.81 E-06 

9.87 E-09 

.04 

0.1398758745 

0.1398758818 

2.80 E-06 

1.57 E-06 

7.27 E-09 

.06 

0.0579829374 

0.0579829414 

1.55 E-06 

8.65 E-07 

4.01 E-09 

.08 

0.0303798097 

0.0303798097 

7.58 E-07 

4.24 E-07 

1.97 E-09 

.10 

0.2274516347 

0.2274516347 

3.49 E-07 

1.95 E-07 

9.04-10 


Table 2: Comparison of Numerical solution of problem 2. 


t 

y(0 

y n +1 

ER 1 

er 2 

Abs 

.2 

1.7755301747 

1.7755301873 

3.37 E-06 

2.13 E-06 

1.25 E-08 

.4 

1.7051273147 

1.7051273208 

1.74 E-06 

1.57 E-06 

6.12 E-09 

.6 

1.6806121251 

1.6806121278 

7.92 E-07 

8.65 E-07 

2.66 E-09 

.8 

1.6717699224 

1.6717699235 

3.44 E-07 

4.24 E-07 

1.12 E-09 

.10 

1.6685404228 

1.6685404233 

1.46 E-07 

1.95 E-07 

4.62 E-10 


Analytic Solution: y(t) — — 
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Table 3: Comparison solution of problem 3. 
Analytic Solution: y(t) — ^ 


t 

y(t) 

y n +1 

ER 1 

er 2 

Abs 

.1 

1.0688853015 

1.0688853144 

9.93 E-08 

5.18 E-06 

1.30 E-08 

.2 

1.6544440817 

1.6544440919 

7.72 E-08 

3.66 E-06 

1.02 E-08 

.3 

2.0180987318 

2.0180987391 

5.55 E-08 

2.49 E-06 

7.31 E-09 

.4 

2.2617841981 

2.2617842034 

4.03 E-08 

1.17 E-06 

5.28 E-09 

.5 

2.4335031421s 

2.4335031460 

2.98 E-08 

1.26 E-06 

3.89 E-09 


IV. Discussion / Conclusion 

The new method has four stages compared with Lie and Norsett with five stages [8] thus the 
implementation cost is cheaper. The new method converges better and faster to the exact solutions (see 
comparison tables above), has biger area of A-stability than methods [3] and [8] consequently the have 
developed a new Gauss- Quadrature Runge-Kutta method with four points or stages only. 
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